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Abstract 



The paper deals with a fractional derivative introduced by means of the Fourier 
transform. The explicit form of the kernel of general derivative operator acting on 
the functions analytic on a curve in complex plane is deduced and the correspon- 
dence with some well known approaches is shown. In particular, it is shown how the 
uniqueness of the operation depends on the derivative order type (integer, rational, 
irrational, complex ) and the number of poles of considered function in the complex 
plane. 



"*-! 1 Introduction 



The fractional differentiation and integration (also called fractional calculus) is a notion 
almost equally old as the ordinary differential and integral calculus. Naturally, each a bit 
more gifted student who has just understood what is the first and second derivative can 
ask the question: well, but what is for example 1.5 - fold derivative? There are several 
ways how to answer such a question. An excellent review on the theory of arbitrary order 
differentiation (generally of complex order) including also interesting historical notes and 
comprehensive list of references to original papers (more than thousand items) is given 
in the recently published monograph |J. Some new results were presented also in the 
recent conference || dedicated to this topic. The fractional calculus has also plenty 
of applications, see also e.g. H],H,[@] and citations therein. Possible use in quantum 
mechanics and the field theory is discussed in |I(J. Recently, fractional derivative was 



mentioned also in || as a particular case of pseudo - differential operators applied in non 
local field theory. 
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Apparently, the general prescription for a definition of fractional derivative is using 
of some representation of ordinary n-fold derivative (primitive function) which can be 
in some natural way interpolated to n-non integer. Actually, following the mentioned 
monograph, all the known approaches are always somehow connected with some of the 
following relations. 

1) The well known formula for n-fold integral 

xi PX2 rx n i «i 

dx 2 dx 3 ... f(t)dt=— {x 1 -t) n - l f{t)dt (1.1) 

Ja Ja r(Ti) J a 

allows the substitution of n by some real a > 0. In this way fractional integration is 
introduced. Then fractional derivatives can be obtained by ordinary differentiation of 
fractional integrals. This is the basis of the construction known as Riemann - Liouville 
fractional calculus. Let us note, that in this approach resulting function even in the case 
of fractional derivatives in general depends on the fixing of the integration limit a on the 
right hand side of (|1.1|) . 

2) The Cauchy formula for analytic functions in some region of the complex plane 

/(•>(„) = 25+1) / «f>* (1.2) 

2tTT J c (z - Zq) n+i 

in principle enables generalization to fractional derivatives, nevertheless the direct exten- 
sion to non-integer values of n leads to difficulties arising from multivaluedness of the 
term (z — Zo) a+1 and the result also depends on the choice of the cut and integration 
curve . 

3) Analytic continuation of derivative (integral) of exponential and power function 

^— exp(cz) = c a exp(cz), ^—(z - cf = ^ + X ) u _ c \p-a ( L3 ) 
dz a PV J PV ; ' dz aK ' T(p-a + l) K J K J 

Obviously these relations allow to define fractional derivatives of the functions which 
can be expressed as linear combinations of power and exponential functions. Also this 
approach is not completely consistent, as can be illustrated by fractional derivative of 
exponential function expanded to the power series 

ex PM = Ef|TT)' ( L4 ) 

but for a— non integer 

d a „ , x , „ (cz) k - a d a /y> (cz" k 

^T(k-a + l) ~ dz a 

k=0 v ' \k=0 



So the task of extrapolation of integer derivative order to arbitrary one has not the unique 
solution. 
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The approach proposed in this paper is based on the fractional derivative of the expo- 
nential function ( |1.3|) entering the Fourier transform of a given function. On this basis in 
Sec. 2 explicit form of the kernel of fractional derivative operator is deduced. In Sec. 3 the 
composition relation for the derivative operator is proved. The generalization of the case 
of the functions on the real axis to the case of function on the complex plane is done in 
Sec. 4, which is concluded by theorem summarizing the results. The last section is devoted 
to the discussion of some consequences following from the theorem and to a comparison 
with the known approaches as well. 

2 Definition of fractional derivative by means of Fourier 
transform 

Let f(x) be a function having Fourier picture f(k) : 

/+oo Z' + OO ^ 

f(x)exp(ikx)dx, f( x ) = 7T~ / f(k)exp(—ikx)dk. (2.1) 

Then let us create the function 

i r+oo ^ 

f<*(x) = — I (-ik) a f{k)exp{-ikx)dk, a > -1 (2.2) 

and define 

i r+oo 

D a ( w ) = — (-ik) a exp(-ikw)dk. (2.3) 
2vr ,/_ 00 

The function f a (x) can be formally expressed 

/+oo 
D a (x-y)f(y)dy. (2.4) 
-oo 



Now let us calculate the integral ( |2.3j ), which depends on the way of passing about the 
singularity k = and choice of the branch and cut orientation of the function k a . For the 
beginning let us assume the cut is given by the half line either (0, — oo) or (0, +oo). For 
complex functions £ a = (£i + i^) a we shall accept the phase convention 



hm (e 1 +^ 2 ) Q =ieri 6,6 >o, (2.5) 
$2^0+ 



i.e. it holds 



£ < £ > cut orientation 

hm (£ + xeY l\r\= I ex P(+ ma ) 1 (°> +0 °) 

mn ^ + ie) 1 1 § 1 < exp(+ ^ a) x (Q) _ Qo) (2.6) 

v ' \a. 1 \ j-a. 1 J exp(+Z7ra) exp(2iira) (0, +00) 

hm (£ - ie) a \ £ a \= < } . { \ } { 

e^o+ y l exp(— ma) 1 (0, —00) 
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Let us define 



D%{w) 



2tt 



hooiiO 



(z/c) Q exp(—ikw)dk. 



(2.7) 



oo±i0 



If we accept phase convention ( |2.6| ) for A; a in the integral ( |2.7| ), then the uncertainty of 
phase of the expression (|2.7| ) is involved only in the factor 



-l) a = exp(ia [2n + l] ir) 



(2.8) 



where n is any integer number. Within arbitrariness given by (|2.8| ) the functions D", Di 
do not depend on the cut orientation and it holds 



D a (w) 








for w < 
for w > 0, 



(2.9) 



which is evident from the corresponding integrals having the paths closed in infinity - 
for D^D ^) in upper (lower) half-plane. Now we shall calculate the integrals ( |2.7| ) in 
remaining regions of w. Let us split them into two parts 



D%{w) 



2tt 



0±i0 



(ik) a exp(—ikw)dk + 



— oo±i0 

and substitute real parameter w by complex one 



hoo±iO 



(ik) a exp(—ikw)dk 



0±i0 



(2.10) 



Z2 



w + le 
w — ie 



for k < 
for /c > 0, 



(2.11) 



where e > 0. This substitution ensures absolute convergence of both integrals in (|2.10| ). 
Next let us make in ( |2.1(J| ) further substitution 



£ = ikzi 
£ = ikz 2 



for k < 
for k > 0. 



(2.12) 



In this way instead of D±(w) we get functions depending also on e : 



:-i) c 

2tu 



rexp(-£K + - 



+1 



^exp(-0^ 



A' 2 



(2.13) 



We assume 

a) Cut of complex function k a is given by half-line (0,+oo). 

b) Function values zf , z% are considered values of one complex function z a in the two 
different points. According to the cut orientation of z a one can get either z% — ( z±)* or 



z")* exp(2z7ra) . Cut of the function z a is assumed 



(0,+oo) 
(0,-oo) 



for E% 
for E a . 



(2.14) 
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integral 


variable 


k < 


fc > 




k 


71 





E% 


z 


(0,7T) 


(7r,27r) 






(37r/2,5vr/2) 


(3vr/2,57r/2) 




k 


7T 


2tt 


El 


z 


(0,7T) 


(-7T.0) 






(37r/2,57r/2) 


(37r/2,57r/2) 



Table 1: The phase correspondence of variables £, z, depending on the cut orientation 



Im 5 




Figure 1: Integration paths in Eqs. ( [2.13j ),( |2.15|) . 
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Latter on we shall come back to these assumptions and judge how they affected our result. 

Assumptions concerning the cuts of k a , z a correspond to phases or to the intervals of 
phases of variables k, z and correspondingly to phases of the variable £. All considered 
possibilities are summarized in Tab.|l| . The corresponding integration paths are shown 
in Fig.|l|. Instead of the interval (37r/2, 57r/2) for arg£ we took the interval (— 7r/2,7r/2) 
since the corresponding phase shift is the same for both integrals in ( |2.13| ) and can be 
included in factor the (— l) a ahead of the integrals. Obviously, it holds 

p poo 

rexp(-£K = -/ rexp(-£K= / f a exp(-f)d£ = T(a + 1). (2.15) 



o 



After inserting into (|2.13|) we get 

„/ s (-i) Q+1 r(«+ 1) f I 1 \ , 

El(w, e = i - V ; — - - — - , e > 0. 2.16 

We shall regard the integrals ( |2.7| ) as the limits 

D a ± (w)= lim£«(«,, f ), (2.17) 



hence the kernel of the operator ( |2.3| ) which can be also considered the generalized function 
is symbolically written 

(-l) a+1 r(a + 1) / 1 1 \ 

where the two modes correspond to different cuts of (w + ie) a+1 

D%(w) for cut (0,+oo) 

Dt(w) for cut (0,-oo). K ' 



Let us note the function ( |2.18| ) is well defined for any complex a = a\ + ^ — 1, — 2, .., 
but for the beginning we assume a2 = 0. Now let us go back to the assumptions a), b) 
which imply result fl2.16j ). 



a) Assuming the opposite cut orientation (0,— oo) for k a in Eq. ( |2.10|) and repeating 
the corresponding sequence of steps does not change anything for E9(w,e) whereas in 
the case of E"(w, e) the phase arg£ will shift by — 2n in both integrals in ( |2.13| ), but this 
change can be included in factor (—1)° ahead of the integral. 

b) Let us assume the cuts in z a having the opposite orientation than in ( |2.14j ). Let us 
take e.g. function E*(w, e), then the phase of £ complies with 

37r/2 < arg£ < 5n/2 for k > , , 

-7I-/2 <arg£< tt/2 for k < 0, ^ ' 



i.e. the phase of second integral in ( 2.13] ) is now shifted by — 2ixa and instead of ( 2.16| ) 
we get 



KM = tS!^£±^(. L_ _ gSfcj^V ,.,„. ,2.2,; 

2m \(w + ie) a+1 (w - ie) a+l J 
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Obviously this function for w G (—00, +00) with assumed cut of z a on (0, —00) is identi- 
cally equal to the function ( |2.16|) with the cut on (0, +00). The analogous result can be 
obtained also for E°_. 

So we have shown, that result ( |2.16|) of the integration (|2.7| ) does not depend on 
the choice of cut orientation of the functions k a and w a . The cut orientation of w a in 
( p,16| ) and ( [2.18| ) is dictated only by the way of passing around the singularity k = in 
initial integrals (|2.7|) and the correspondence (|2.19|) always holds. Next let us notice the 
important property of the functions D±, E±. Obviously it holds 

d -E£(w,e) = El + \w,e), 



dw (2.22) 
(> Dl(w) = Dl +1 (w). 



dw 



Now let us go back to the Eq. (|2.4j) and consider how to calculate this integral. It is 
possible either 

a) First to make the difference of both terms in (|2.16|) , then integration and limit for 
e - 0, 

or 

b) First to calculate both integrals independently, then make the limit of their differ- 
ence. 

Let us discuss the both ways and compare the results. 

a) The action of operator D± as the integral of difference 
We make calculation separately for the two cases: 

al) a = n > is an integer number 
In that case the complex function w a has no cuts (we can omit the subscript ±) and for 
n = 0we can write 

E\w,e)= l --^— 2} (2.23) 
i.e. for e — > we get the known representation of 5-function (see e.g. J3|,p.35) 

5(w) = lim - ■ (2.24) 

which acts on the function /, 

/+00 1 r+00 

^ S(x - y)f(y)dy = hm - J ^ j- - y? - - f(y)dy = f(x). (2.25) 

Using equations (|2.22|) — (|2.25| ) one can easily show 

= lim / E n (x- y, e)f(y)dy (2.26) 

dx n e^0+ t /_ OG 



7 



for any n > for which the integral converges. So it is possible to identify 



D n {w) -- 

d n f(x) 
dx n 



= J^D n (x-y)f(y)dy, 



i.e. action of the operator D n corresponds to n-fold derivative. 

a.2) a is a real, non integer number 
Taking into account the cut orientations (|2.19|) calculation of limits (|2.17 ) gives 



D!(«j) 
D+(w) 

D a (w) 








T(a + 1) sin([« + 1]tt) 



,a+l 



D%{w) = + 



TTW 

T(a + 1) sin([a + l]vr) 



TTW 



a+l 



w > 

w < 

w < 

u> > 0. 



(2.27) 



(2.28) 



The first two equations consist with Eqs. ( |2.9| ), therefore again confirm correct correspon- 
dence of cuts in ( p.!9|) . After inserting of ( [2.28j ) into Eq. ( p.4|) we get 



F(a + 1) Bin([a + l]vr) 



7T 



=poo 



f(y)dy 
(x-y) 



a+l 



(2.29) 



For y = x and a > the integral has a singularity. The method for regularization of this 
integral will become apparent in the next. 

b) The action of operator as a difference of two integrals 
After inserting of Q2.18D to ( J2.4Q we get 



lim 



r(« + r 



2i7T 



f{y)dy 



(x — y + ie) a 



+i 



f(y)dy 



[x — y — it 



\a+l 



(2.30) 



We assume for the present the function f(y) is analytic on the whole real axis. Last 
equation can be rewritten 



flix) = lim 

J±y ' 6^0+ 



\a+l 



r(« + i) 



2in 



f(z)dz 



\x — z 



I a+l 



f{z)dz 
(x — z) 



a+l 



(2.31) 



Again let us separate two cases: 

61) a = n > is an integer number 
In Eq. ( |2.31|) we can link up both integration paths in z = ±oo and write 



lim 

e^0+ 



l) n+1 r(n+l) f f(z)dz d n f 



2m 



c 



(x - z) n+1 



dx n 



(2.32) 
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Im z 



Li- 



Re z 



Figure 2: Integration paths in Eq. ( p.33j ): a)C + , 6)C_. 



where C is any closed curve enclosing the singular point. Therefore for n > the operator 
D n can be again identified with ordinary n-fold derivative. 

62) a is a real, non integer number 
Similarly, as in the case 61) we get the integrals 



/±(*) 



\a+l 



lim 



r(« + i) 



2ivr 



f{z)dz 



c± 



(x — z 



(2.33) 



where the paths C± can pass about the corresponding cut as shown in Fig|| The integrals 
converge provided that 



lim M 



0. 



(2.34) 



The integration paths can be decomposed into three parts L 1± , L 2 ±, L 3 ± and after eval- 
uation of the corresponding integrals one gets 



lim 



r(a + 1) sin([a + 



71 



f{z)dz 
(x — z) 



a+1 



+ 



(2.35) 



or using the relation ( [A.lOp from Appendix, 



lim 



-o+ T(-a) 



f(z)dz , (±l) a f(x) 



(x — z) 



a+1 



+ 



(2.36) 



For a < the integrals are finite and the second terms vanish. In this case Eq. ( 2.35Q is 
identical with Eq. ( |2.29| ). On the other hand we have assumed the initial integral ( |2.33| ) 
is finite, therefore the sum of both terms in (|2.35 ) is also finite even for any a > 0. In this 
sense by addition of the second term the integral in Q2.29| ) can be regularized. Note, that 
in contradistinction to ( 2.18 ) the relation ( |2.36| ) is well defined also for a = — 1, —2, —3, ... 
but a ^ 0,1, 2, 3, ... The last equation for a negative integer 



Tin) 



[x — z) 



\n-l 



f(z)dz 



(2.37) 



+ OO 
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is a special case of the n-fold integral formula fll.lj) . 

Now, let us assume < a < 1 and calculate integral (|2.33|) by parts. Obviously 

=poo± iO 



f(z)dz 
(x - z) a+l 



1 

a 



f'(z)dz 



+ 



(x — z) 



(2.381 



=poo=F iO 



'c± {? ~ z ) 

and integral on right side is finite. Any a can be written as the sum of integer and 
fractional part 

a = n + Aa, n = [a], < Aa < 1. (2.39) 

For n > we can repeat integration by parts n + 1 times and if the function / meets 
requirements 

f p (z) 



lim 







»±oo (x - z) a ~P 

then instead of (|2.35| ) we get 

[-l) Aa T(Aa 



lim 

6^0+ 



for p = 0,1, ...n, 
f n+ \z)dz 



(2.40) 



2m 

r(Aa) sin(Aa7r) 

7T 



x — z 



lAa 



T-s. fn+l^ dz 

(x — z) 



Aa 



(2.41) 



Let us note that integrals ( |2.33| ) can be modified also in other way. We assume that all 
integrals 

f(z)dz 



I±(l,x)= [ 

JC- 



, c± (x-zp' 
converge, then recurrent relation holds 

d 



7 = Aa, Aa + 1, ...a + 1 



dx 



(2.42) 



(2.43) 



Application of this relation for integral ( p.33| ) and a = n + Aa gives 

_ \Aan/ a _ \ / j \ "+ 1 

= lim ' 



-l) Aa r(Aa) 



2Z7T 



dx 



f {z)dz 



(x - z 



Aa 



r(Aa) sin(AaTr) ( d \ n+1 f T °° f (z)dz 
dx) 
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(x - z 



Aa ' 



(2.44) 



Actually in Eq. (|2.41 ) we apply operation D Aa_1 on n+1 - fold derivative of the function 
/ whereas in Eq. ( |2.44| ) the both operations are interchanged. Using relation ( A.10|) 
obviously we can modify both equations: 



rri- Aai 



d 



f n+1 {z)dz _ 1 
;poo (x - z) Aa ~ r(l - Aa) \dx 



n+l 



f {z)dz 



=poo 



(X — Z 



,Aa ' 



(2.45) 



So we have shown the operator D Q with the kernel ( |2.1<j| ) can be considered a contin- 
uous interpolation of the ordinary n-fold derivative (integral) of the functions analytic on 
real axis and fulfilling the condition ( 2.34 ). 



10 



3 Composition of fractional derivatives 



In this section we shall investigate how the composition of the operators D a is realized in 
the representation given by Eq. ( j2.18|) . Therefore we shall deal with the integrals 



-oo 



D a {x-y)D^y-z)dy. (3.1) 



Let us denote 



(W + \ T)7 r>0 (3.2) 



(w — ir)' 



and 



h = h'(a + l,x- y)h\(3 + l,y- z)dy 
h = f^™ h,(a+l,x- y)h.(f3 + l,y- z)dy 
h = f^™ h.(a+l,x- y)h'(f3 + l,y- z)dy 
h = f-™ h'(a+l,x- y)h.(f3 + l,y- z)dy, 



(3.3) 



then 



4tT 2 

Now let us calculate the more general integral 



-iW^ + ijiv + i) > + /j + 



J= I 77 ^Jl 7TS7 . (3-5) 



K 



(z 2 - z ) a + i (z- Zl y+ r 



where K is arbitrary line in complex plane, Zi, z 2 any two (diverse) points and let a + f3 > 
— 1. We also assume the cuts of the function in the integral do not intersect the line K. 
Then there are two possibilities: 

a) K is not passing between the points zi,z 2 , see Fig.^|a. Then obviously J = 0, 
since the line K can be closed in infinity by the arc in half plane which do not contain 
singularities z±, z 2 . 

b) K is passing between the points z±,z 2 , see Fig.^. K' denotes line crossing the 
segment (zi, z 2 ) perpendicularly at its center. If we assume, that the cuts do not intersect 
any of both lines K, K', then in the integral ( |3.5| ) path K can be substituted by K' . 
Further, if we denote 

zq = (zi + z 2 )/2, rexp(z^) = (z 2 - zt)/2 (3.6) 

and substitute z = z + itexp(iip), then 

i C +0 ° dt 

ex.p[icp(a + (3 + 1)] 7_ 0O (r - it) a+1 (r + it)P +1 ' v ' ; 
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Figure 3: Integration paths in Eq. (|3.5|): a)case when the integral vanishes, 6)case leading 
to the result ( gjg) . 



The last integral can be found in tables (see e.g. 0,p.3Ol), using (|3.6| ) we obtain 

7 _ 2Z7T r(a + /? + l) 

(z 2 -z 1 )^+ 1 T(a+l)T((3 + iy 1 ' ) 

Let us note that the opposite orientation of the line K (i.e. point z 2 on left side in respect 
to the direction of K, as the line L in Fig.|3|&) should give result ( j3.8|) with opposite sign. 
Obviously the integrals I3, J4 vanish and for I 1; J 2 we get 

= -2m T(a + P+l) 

1 (x-z + 2zr) Q +^+ 1 r(a + l)r(/5+l) 1 ' ' 



+2i7r r(q + /3+ 1) 

2 ~ (a; - z - 2ir)«+/ 3 + 1 r(a + l)T(/3 + 1) 

and inserting into (|3.4j ) gives 



(3.10) 



j_ (-i) Q+/?+1 r(q + p + f 1 1 \ f3ir 

2ivr \(x- z + 2ir)°+^+ 1 (z - z - 2ir) a +^ +1 / ' 1 ' ' 



For r — > this expression corresponds to D a+l3 (x — z) in ( [2.181) . 

This result is formally correct, nevertheless its drawback is that it does not reflect the 
correspondence of cut orientations in initial expressions ( |3.1D and the final ( ^.11[ ). More 



rigorous discussion about the cuts we postpone to Appendix, here give only result: the 
following composition relation holds for operators ( [2.18|) with equally oriented cuts 

/+00 
D%(x - y)Dl{y - z)dy a, /3 ^ -1,-2,-3...; a + /3 > -1. 
-00 

(3.12) 
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Let us note, that validity of the relation ( p-12j) can be verified in the initial representation 
(]2. 71) as well. Further, considering D± as a generalized function, then the condition 
a + (3 > — 1 can be omitted and the composition relation has the form 

/+oo r+oo r+oo 

D a ± + P(x-0f(0dZ= / D%(x-y)Di(y-t)f(Z)dtdy, a, {3 ? -1, -2, .. 

-oo J — oo J —oo 

(3.13) 

for all functions analytic on the real axis and fulfilling 

lim ^S = 0. (3.14) 

Eq. ( |3.13| ) follows from ( |3 . 1 2| ) and relation -^-D'jT 1 (w) = D±(w). Repetitional integration 
by parts on both sides can reduce sum a + (3 by any natural number, so in this way 
validity of ( |3.13| ) is proved. 

All our previous considerations concerned acting of the operator D a on real axis. In 
the next we shall try to enlarge them on whole complex plane. 



4 Fractional derivative in complex plane 

First, let us illustrate the notion fractional derivative introduced in previous part on a 
particular example. Take the function 

/M = —> = --(—-—), (4-1) 
1 + x l 2i \x + i x — i J 

for which (|2.33|) gives 

\ r (-l) a+1 I> + l) f dz fA , 

fl(x) = hm / — — --. (4.2) 

J±K ' e-^o+ 2m J c± (x-z) a+1 (l + z 2 ) K ' 

Obviously for a > —2, the integration paths in Fig.|2] closed by arcs in infinity as shown 
in Fig.f|, are possible and 

f « (x) (-ir +1 r(^ + i) / i i \ 

Even though the result formally does not depend on given subscript (+ or — ), it is 
necessary to take into account different cut orientation for /?, In accordance with the 



phase convention i 


m 


we take 


(x + i] 


a+l 


R a+1 


(x-f 


a+1 


R a+1 


(x + i] 


a+l 


R a+1 


(x - l 


a+l 


R a+1 



exp[i(27r — <f)(a + 1)] 



for/- 

(4.4) 

for /« 
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Figure 4: Integration paths in Eq. (|4.2|) closed in infinity: a)C + , 6)C_. 

where 



R = vl + a?, 99 = — — arcsin — . (4.5) 

2 R 



After inserting into (|4.3|) and a simple rearrangement, we obtain 

a+l 



ivrh) T{a+1) sin [ ( arcsin w+w ± l) (a + 1] 



(4.6) 



Remark: From the last relation it can be easily shown that e.g. 

a) for a integer, a = n > — 1, it holds /?(#) = f1(x) = f n {x), where / n is ordinary 
n-fold derivative 

6) = (-l) Q /^(-x) 

c) for a — > — 1 we get a primitive function to / : f±{%) —> arctan(x) ± 7r/2. The 
same results can be obtained also from (|2.45|) for n = —1, Aa = 0. 



Using formula ( |2.33| ) we shall try to generalize the operator of fractional derivative on 
real axis to the whole complex plane. For this operator we shall demand again: 

V a f(z) = -± for a = n > 0, (4.7) 



D tt o = B a+lS . (4.8 

For making the generalization more transparent, we shall do it in several steps. 

1) Let us go back to Fig|| and put a question what will change in results ( |OD — (|4.6|) 



if the cut (and integration path correspondingly) would be oriented otherwise than along 
real axis, but e.g. as shown in Fig.^. Obviously for Eq. (|4.3|) nothing will change, but 
the form of Eqs. ( |Q1) , (|4.6D will depend on the mutual position of the cut and both poles. 
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Figure 5: Tentative modification of integration paths in Eq. (|4.2|) and Fig 
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Figure 6: a)Integration paths in Eq. ( [2.33| ) for the function (fO|). (fi, ^2 are phases 
corresponding to terms w% +1 = (zq — Zk) a+1 in the result of integration ( |4.11| ). b) The 
same phases represented for variable w. 



Because we have accepted phase convention ( |2.6| ) only for cuts on real axis, it is now 
necessary to make a more general consideration to determine the phases of both terms in 
Instead of (J4.1|) let us take a complex function 



/(*) 



a 1 



+ 



a 2 



Z — Z\ Z — Z2 

and integration paths C displayed in Fig||a. The corresponding integral will be 

a-y CL2 



/"(*>) 



-i) Q+1 r(a + r 



(z - z,)^ 1 + ~(z - z 2 ) Q+1 



(4.9) 



(4.10) 



The path C passes about the cut of function 1/w 



a+l 



1/(zq — z) a+ in variable z, i.e. 



cut of the function l/w a+1 is oriented in opposite direction, see Fig.^. The phases of 

must be fixed in respect to this cut. Fig.||& 



a+l 
Wi 



z - z x ) a+1 and < +1 
prompts the following rule. 



1 a+l 
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Rule 1: Phase of complex variable w is given by the angle ip of arc leading from positive 
real half axis and measured in positive direction (against clockwise sense) to the point w 
and if the arc intersects the cut, ip is reduced by 2n. 

This rule can be applied also directly for situation in Fig||a for fixing phase of (zq — z) a+1 . 
The only modification is that (p is measured from half line (zo, zq — oo). Therefore angles 
(px, <P2 m Fig.|| fix phases in ( [4.1 0|) 

f (z ) = (-1) F(a + 1) ^ |(, _^ + i| + \(z -z 2 )^\ ) ■ (411) 

Now the introduced rule can be applied also for integration on paths Ci(C 2 ) in Fig.[5| 
Doing this, we shall get the same result as in the case of integration on paths C + (C_) in 
Fig.^. Therefore depending on the position of chosen derivative cut in respect to the both 
poles, the function ( has (up to factor fl2.8| )) in a given point x two different values 
of fractional derivative given by Eq. (|4.6| ). 

2)We shall now generalize the prescription for calculating of fractional derivative of 
function having the form (|4.9|) for functions having a finite number of poles 

TV 

ft w = E (g -V- ni ~°- (412) 



Let the derivative cut be given by half line L = [zo, zq + exp(z#)oo] which does not go 
through any of poles z k . Then 



ton (-D^ + Df /" 

y °> e ^o+ 2m J C{L) {zo - zY+\z - z k )^ 



N 

E 

fc=i 



T(a + n k + l) a k exp(-iip k [a + 1}) 
r(n fc + l) \z -z k \ a+1 (z -z k ) n * 



(4.13) 



integration path C(L) is shown in Fig.0. Angles ip k are calculated using the Rule 1, 
therefore it is obvious, that the function h a (zo) in (|4.13|) can have, according to the cut 
orientation, as much values, as much different poles z k the function (|4.12 ) has. 

3) Now let us consider functions (4.12), p.l3|) in the case, when the cut is not a line, but 



some general curve connecting the points zq and exp(i#)oo. Then the result (|4.13|) will be 
formally the same, only the prescription for fixing of angles <p k must be modified. General 
cases are illustrated in Fig.||a,6, where it is apparent that the arc, which 'measures' angle, 
may intersect the cut several times. The consistent assigning of angles ip k corresponding 
to points z k in Eq. ( [4.13| ) can be ensured by the following prescription. 
Rule 2: Let us choose on the half line (zq, zo — oo) any reference point zr. Angle <p k 
is given by angle ZRZoz k measured in positive sense and then for each intersection with 
the cut is corrected as follows. Superpose palm of right or left hand at an intersection in 
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Figure 7: Integration paths in Eq. ( |2.33| ) with the function ( [4.12| ) and phases ifk appearing 
in the result of integration (}4.13|) . 
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Figure 8: a), b) Examples of the curvilinear cuts generating integration paths for Eq. 
Q2.33| ) with the function Q4.12 ). Other symbols are defined in Rule 2. 
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such a way, that fingers lead from zr towards z k and thumb leads in direction of the cut 
from branching point z . If this condition is met by right (left) hand, tpk will be enhanced 
(reduced) by 2n. 

Remark: It is substantial for all Zk to choose one common reference point zr. A shift 
of this point results at most in only equal shifts of all angles <fk — > y?fc + 2mr, which do 
not change Eq. ( |4.13| ). 



Therefore for any curvilinear cut the relation ( [4.13 ) can be written 



h * (z ) = ( 1 x a y> r(a + 7i fc + l)a fc exp(-i[y fe + 2m fc 7r][a + l]) 

• -° j ~ { ' ^ T(n k + 1) \za-z k \ a+1 (z - z k )^ ' ' 



k=i 



The set of integer numbers (or more exactly their differences) characterizes the way how 
the curvilinear cut passes among the poles z k - If we accept curvilinear cuts for derivative 
of the function ( |4.1| ), then using ( |4.14| ) one can obtain 



fa, x (-l) fc(a+1) r(a + l) . 
f (x) = — — —— — Ti — sm 

(VT+x^) a+1 



( arcsin 7x=f + ( 2fc+1 )i) 



(4.15) 



where k = ni2 — m\. Let us note that 



lim f a (x) = arctan(x) + (2k + 1)-, (4.16) 
a— »— 1 2 



i.e. we obtain the infinite (but countable) set of primitive functions for the function (4~T). 



4) So far we have considered only analytic functions of the form Q4.12 ), for which the 



corresponding integral on whole circle C(oo) vanishes. Now we are going to consider the 
general case, when this integral vanishes only on a part of the circle. But first, let us go 
back to the operator ( |2.18| ) acting on analytic functions on real axis (or the part of this 
axis) and try to generalize it for analytic functions on a curve in complex plane connecting 
a pair of points on the circle C(oo), see Fig.|9|a. Let this curve be given as a continuous 
complex function of real parameter t G (— oo, +oo) 

z = x(t) + iy(t) = ip(t), ip(—oo) = ooi, ^(+oo) = oo 2 , (4-17) 

then its derivative 

determines in complex plane the vector tangent to curve ip and oriented in direction of 
increasing t. The function 

u(t) = 7=ff2=^ = expM*)] (4.19) 
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Figure 9: a)Integration paths for operator ( |4.20| ) in complex plane, b) Corresponding 
curvilinear cuts in Eq. (|4.20|) . c) Corresponding integration path in Eq. (|4.23|) . 



represents this vector in its normalized value (uj(t) is phase of this vector) and the function 
iv(t) normalized vector perpendicular to tp and oriented left in respect to the course of i/j. 

Now let us define the function of the complex variable z6i|i with complex parameters 
z G tj) and a ^ —1, —2, —3, ... 

-l) a+1 T(a + 1) / 1_ 

♦o+ 2m \ (z — z + ieu) a+l (z — z - ieu) a+1 

having the curvilinear cuts (for a ^ 0, 1, 2, ..) corresponding to the both terms coming out 
the points z§±iev and going jointly along the curve ip to points ooi or oo 2 , see Fig.^. So, 
in contradistinction to linear cuts the form of the cut of w a+1 = (z Q — z) a+1 does depend 
on position of point z on ip. 

Next, using this function we define the operator 



D°L(zo - z) = lim ^— ^ 1 M r^r (4-20) 



/ = /; ± (z ) = / B%(zo - z)f(z)dz = / D^(t ) - j>{t)) f m)W Wt (4.21 



acting on the functions analytic on the curve ip for which the integral converges. 
Remark: Let us note the integral (|4.21| ) depends on choice of cut end-point (ooi or oo 2 ) but 



does not depend on in which direction (ooi — > oo 2 or ooi <— oo 2 ) integration is done. That 
is a consequence of the fact that change of integration direction ifj(t) — > tp(—t), dz — > —dz 
implies the change v{t) — > —u(t), which implies the change -D^(z — z) — > —D^(zo — z), 
therefore the integral does not change. 

If we calculate ( |4.21|) as the difference 
D7 f = lim + !) ([ fWz f M± ^ (4 22) 

then this difference can be expressed as 

r D«(z z) f(z)dz - lim i-^^ + l) [ f(z)dz 

I D4z ° - z)fiz)dz - i& — ^ — y c( ,, e) (,o-^ ! (4 - 23) 
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Figure 10: Domain G (grey area) and the curves defined in assumptions of the Theorem. 



where the path in 'distance' e passes about the cut coming out branching point z on the 
curve if; to infinity, see Fig.||c. 

To label the integration cuts ending either at ooi or 002 we can accept the following 
convention. Let 00 1 = exp(z6 l i)oo and 002 = exp(i6 , 2 )oo, then for 



„ , n \ COS 6^ < cos 6 2 
COS 61 7^ COS 62 < Q ^ Q 

[ COS 0i > COS 6*2 

„ f sin 61 < sin 9 2 

COS 61 = COS U 2 { . n . n 

\ sin X > sin 9 2 






= (z ,ooi), 




= (Zq, 


oo 2 




= (z ,oo 2 ), 


i>- 


= fa, 


00 1 




= (zo, OOi), 




= (zq, 


oo 2 




= (^o, oo 2 ), 




= {zq, 


00 1 



(4.24) 



and correspondingly we index related symbols D^ ± , D^ ± , f%±, C±(if),e). 

Let us note the definition of the operator by Eqs. (|4.20|) , (|4.21|) ensure that the 
both curves C±(if), e) are oriented in such a way that after their closing by the circle C(oo) 
see Fig.|]c, there arises closed curve having always clockwise orientation. We denote these 
curves as K±(ip, e). 

Now, after all these preparing steps we can come to the formulation and proof of the 
following theorem. 
Theorem 

Assumptions: 

i) G is the domain in complex plane containing the part (or parts) C\ = Gfl C(oo) of 
the circle C(oo). The curve C\ forms one part, the curve Cq is the second part and closed 



curve Cq = C U C\ constitutes the complete domain boundary (Fig. |lCi|) . 
ii) if) is a curve defined in (|4.17f) and if) C G, if) fl C = 0. 

Hi) Function f(z) is analytic in G except for the poles z^ £ if), z^ ^ Cq, k = 1...N 
and for a given a = a± + ia 2 7^ —1, —2, —3... and any z\ G C\ meets the condition 

f(z) 

lim = 0. (4.25) 
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Statements: 
1. Operation 



D£t/ = f^(z ) = [ D« ± (z - z)f(z)dz, z,e^ (4.26) 

Jib 



is a fractional derivative satisfying the conditions ( |4.7|) , (|4.8| ). In particular that means 
the operation does not depend on if) if a = n > is integer. 
2. /£ ± is given also equivalently by 

l) a+1 r(a + l) f f(z)dz 



f$±(zo)= lim v ' v y / JV y (4.27) 

3. For iV = the value f2±{za) does not depend on if), i.e. up to factor (—1)° is 
uniquely defined. 

4. For N > 1 and a non integer there remain cases: 

4.1. «i is rational and a 2 = 0, then f% ± (z ) does depend on ■0, but the set of values 
is finite. 

4.2. a.\ is irrational or a 2 ^ 0, then f$±(zo) does depend on if) as well and in general 
the set of values is infinite, but countable. 

Proof: 

We shall start from the statement 2. Obviously, its validity follows from Eqs. (|4.22|) , 

US). 



Now let us consider the statement 4. If we express the function f(z) as the sum 
f(z) = g(z) + h(z), where g(z) is analytic in G and h(z) has the form (|4.12|) , then 
obviously 

r (z)- hm i-^^ + i) ! f^ dz 

= l im (~l) a+1 r(« + l) / f(z)dz _ (-l) a+1 r(q + l) /■ f(z)dz 

~°+ 2Z7T A' ± (V>,e) Oo - Z) a+1 2Z7T ^ (zq " z)^ 1 

-i -wa. I r ( a + nfc + X ) a fcexp(-i[y fc + 2m k 7i}[a + 1]) | 

Vt^ r ^ + 1 ) I *b " l a+1 (*b " + ^ 

where we have denoted 

Hzo) = r( " + 1) / ? +1 • (4.29) 

The angles + 2m k 7i corresponding to the poles z k are evaluated according to Rule 2. 
In the function I(zq) the orientation of the path Cq is accordant with K± and phase of 
(20 — z) a+1 must be also according to the Rule 2 related to the same reference point zr. 
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Now, if «i = p/q (p, q are not commensurable) and «2 = 0, then each term in the last sum 
has the same value also for m' k = + q, i.e. depending on the form of the cut and the 
corresponding set {m^, 1 < k < N} the expression ( 4.28 ) has only finite number of values 



for some zq. On the other hand for a irrational different sets {m^} give different values 
of the sum, therefore in general the number of values is infinite. For «2 7^ multi-value 



factors in the sum ( |4.28| ) are expanded 

l a = exp(2z7rmo;i) ■ exp(— 2nma2) (4.30) 

i.e. mai,m«2 determine phase and scale of individual terms in the sum. Obviously for a 
complex a the set of values f a {zo) depending on ip is in general infinite, similarly as for 
a irrational. Therefore the statements 4.1, 4.2 are proved. 



For N = Eq. (H~28|) is simplified: 

f$ ± (z ) = (-l) a I(z ) (4.31) 

and trueness of statement 3. is evident. 

Finally, let us consider the statement 1. Validity of the condition (f4.7|) follows from 
the statement 2. For a integer the path C±(ip, e) can be closed around the pole z (having 
positive orientation) and we get the Cauchy integral 

H] 1 f{z)dZ -r(zo). (4.32) 



2m J c(zo) {z - z )" +1 
The condition ( |Q| ) requires validity 

/ D$f(x-Z)f(Z)dt= [ [ D^ ± (x-y)D^ ± (y-Of(0^dy. (4.33) 



This relation follows from ( |4.21|) and (p,13|) in which the substitutions £ — > ^(0 an d 



y — > ip(y) are applied. So the whole proof is completed. 

5 Discussion 

5.1 Some remarks regarding the theorem 

Now, let us look into the statements of the Theorem more comprehensively. The function 
I(zq) in ( |4.28| ) can be expressed 



r(a+l) f f{z)dz 



2m J Co {z - z)«+ 1 
r(a + 1) exp(— i2mii[a + 1]) f f(z) exp[— iip(z)(a + l)]dz 



2™ J Co l(^0-^) a+1 | 



(5.1) 
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Since the cut if> does not intersect the curve Co the phases of all points on Co are 'corrected' 
by the same factor standing ahead of the integral. Now it is obvious that for a = p/q 
the function (|4.28| ) can have at most q N+1 different values, including multi-value factor 

Actually, the function (|4.28|) can be written as 



f(! , , _. u n I \^ r(o + 7ifc + l) a k r(q + l) f f(Qd£ ] , 

/lJ ' 1 J ^ r(n* + i) (z-^+«+ 1 + 2m y Co (* - O"* 1 ' 



and considered a multi-value function with the value at point z depending on the choice 
of cut if> in ( |4.26j ). At given point z there are equivalent any two cuts for which the 
region closed by them and the curve C\ does not contain any pole Obviously, except 
for the points z& the function f a (z) is (as original function /) analytic in the domain 
G. Moreover, in the region of a in which f a = ~D a f exists, this function is apparently 
analytic also in respect to a. 

Special case takes place when a is a negative integer. Then due to the singularity 
r(— n) the kernel (|4.20| ) loses the sense. Nevertheless, assuming in ( |4.20|) initially a 7^ —n 



one can proceed to representation ( [4.27 ), then making the limit e — > (like Eq. (|2.35|) ) 
gives 

r(a + 1) sin ([a + f f f(z + VQe)dz f( z o) 



fS±{zo) = - lim v u J ' \ " ' . + "r^V , (5-3) 

where ^0 is given by Eq. ( |4. and represents direction of the cut at zq (orientation is 
assumed from zq to infinity). This formula already makes sense for a = —n (if condition 
fl4.25|) holds). Using relation (|A.10|) gives 



/£(*b) = / («> - z) n - l f{z)dz, (5.4) 



which is a modification of the known formula (|1.1| ) for n-fold primitive function. The set 
of values f~ n {zo) depends on if) as follows. If f(z) has no poles in the domain G, then 
any two integration paths if>i , 1P2 in (|5.4|) can be connected by a fragment of C\ , in this 
way there arises closed curve C and it holds 



i.e. / n (z ) is determined uniquely. Now, let us suppose f(z) has inside curve C just one 
pole, for z — > z p 

f{z) - (5.6) 

w (z - z p ) n p +1 v ; 



Then obviously any curve C in the integral (|5.5|) can be always substituted by a couple of 
curves Kq, each of them is closed in the phase range (0, 2ir( and having the pole z v inside 
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(e.g. circles centered at z p ). Instead of ( |5.5|) one get the difference 



(n - 1)! 7 K) (2 - z„)" r +' 

where m is the integer depending on the shape of the curve C (m represents the number 
of 'twists' on the C). Using the Cauchy formula, the last equation gives 

n p > n 

A p ={ 2iirma p (5.8) 

{n-n p -l)\n p \ ( *° Zp) "* <n 



For more poles the all corresponding terms (|5.8j) are simply added. For example integra- 



tion constants in (|4.16|) representing the case n — 1, n p = fulfill ( |5.8|) . The composition 



relation for a = — n, /? < and ^ =real axis is proved in Appendix (Eq. (|A.17|) ) and 
apparently can be transformed to any curvilinear path if). 

Let us note, for a negative integer only the representation of D" given by Eq. (|5.3|) (or 
( |5.4| )) makes sense and conversely for a non negative integer only the representation given 
by Eq. ([4.26 ) with kernel ( 4.20Q (or equivalently by Eq. ( 4.27Q ) is well defined. For any a 



complex but non integer both representations are well defined. These two representations 
differ in the corresponding integration paths: 

i) Eq. ( |4.27| ) - integration on some curve enveloping the cut, the path can be closed. 

ii) Eq. (|5.3| ) - integration on the cut itself, the path cannot be closed. 
For a integer it is specific, that the cuts disappear. 

Perhaps most restrictive assumption in the Theorem is the condition ( |4.25|) . The 
question, if one can in some consistent way apply D a to the functions not obeying this 
condition in any part of C(oo) (and simultaneously having the ordinary derivatives - 
primitive functions) requires a further study. Obviously one possible way is to consider 
such functions as generalized functions as well. 

5.2 Concluding remarks 

Have we said something new? First let us show what is not new. Apparently: 

a) The content of Eq. ( |2.45| ) is almost identical with Liouville definition of right and 
left -handed fractional derivatives in ||, p. 95. The only distinction is in phase of /" since 
in the cited definition for the real functions the real value of the derivatives is ad hoc 
postulated. 

b) Eq. (|4.27|) is the Cauchy type integral which ordinarily serves as one of possible 
starting points for the definition of the fractional derivative in complex plane, see ||,p.415. 
In our approach the integration path is uniquely defined by chosen cut. 

On the other hand, the new seems be the following: 

1) The general form of the kernel ( f4.20|) from which the both above mentioned formulae 
follow. 
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2) The construction based on the integration paths enveloping curvilinear cuts, which 
in the result allows to identify fractional derivative - integral with the multi-valued func- 
tion and to determine how the number of values depends on the derivative order type and 
the number of poles which the given function has in the considered region. 
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Appendix: The correspondence of cuts in the composition relation 

Similarly as in (|3.2|) we define 



(tt ' + ; T)7 r>0, (A.1) 



- in-Y< 



[W — IT 



where indices ± denote the cut orientations (0,±oo). Let us consider the integrals 

/+oo 
h(a, x — y)h(b,y — z)dy (A. 2) 

-00 

for various combinations of the cut orientations and their locations above or below the 
real axis. First let us assume 

a < 1, b < 1, a + b > 1. (A.3) 



For the more general integral ( |3.5| ) we have shown the integral vanishes, when the in- 
tegration line does not separate off the both singularities, which is the case of the eight 
integrals ( A.2p involving combinations h ± h±, h^h^, h±h ± , h T h^ . Let us calculate (A.2) 
when e.g. x < z and the both singularities are above the real axis. Then 



|/| = |exp(m7r)/i + I 2 + exp(-ibTr)I 3 \ = 0, (A.4) 



where 



dy 

1 1, = I _ _ L x = (-oo,x), L 2 = (x,z), L 3 = (z,+oo). (A.5) 



Using simple substitutions in the known relation 



x 



A-l, 



r(A)r( M ) 



r(A + ^) 

and denoting d a+b ~ l = \(x — 2;) a+6_1 | one can get 

1 r(l-a)r(a + 6-l) 



1 - xY^dx = ) J VP . ; (A.6) 



h 



d a + b ~ l T(b) 



i r(i-o)r(i-6) 



d a+b-i r(2 - a - b) 

1 r(l - b)T(a + b - 1) 
d a+b ~ l Via) 



(A.7) 
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After inserting into ( |A.4j ) we get 

r(l-a)r(a + 6-l) r(l-a)r(l-6) 



cos (air) 



T(b) 



+ 



r(2-a-6) 



, r(i-6)r( a + 6-i; 

+ COS(07T) = 



r(a) 



(A. 



sin(a7r)r(l — a) sin(&7r)r(l — 6) 



(A.9) 



r(6) r(a) 

The last identity also follows from the known formula (see e.g. [Q], p. 256) 

T(7)r(l -7)sin(7vr) = 7T (A.10) 



Now let us calculate ( A.2j) for remaining combinations ft ft*, h^h^, h±h±, h T h±. For 
r — > it holds 



lim ft, (7, ty) 



expu7rc) 



(A.11) 



where the phases c of functions h entering the integral (|A.2 ) are in accordance with the 
convention (|2.6| ) summarized in Tab.||] This table enables to obtain the phases Ck of their 
products which are listed in the first three columns of Tab.[| The integrals of the all 
combinations summarized in the table can be expressed like ( [A.4Q , 



I = y^exp(ic fc 7r)/ fe 



fc=i 



If we denote 



G = 



2m T(a + b-l) 



d a+b - 1 T(a)T(b) 



(A.12) 



(A.13) 



then using identities (|A.8| ),( |A~T0| ) the sum (|A.12|) can be evaluated. The results are given 
in the last column of Tab.[3| Let us compare the corresponding rows in upper and lower 
part of the table. Obviously for any x, z one can write 



lim / l^fax -y)h ± (b,y - z)dy = lim — ft ± (a + b - 1, x - z) 



and equally for h±. So far we have assumed (|A.3|) , however using identity 

d 



dw 



ft (7, w) = -7/1(7 + 1, w) 



we can enlarge validity of (|A.14 ) to any a, b 

a + b>l, 0,6=^0,-1,-2,-3,... 



(A. 14) 



(A.15) 



(A.16) 
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h(a, x 


-y) 


h(b, x 


-z) 




y < x 


y > x 


y < z 


y> z 


h + 





— a 


-b 





hr 





— a 


-b 





K 


-2a 


— a 


-b 


-2b 


h_ 





a 


b 






Table 2: The phases c depending on the cut location 





y < x 


x < y < z 


z <y 


I 


h + h+ 


-b 


—a — b 


—a 


+ exp(— m\a + b])G 


h + h + 


-2a -b 


—a — b 


-a - 2b 


— exp(— m\a + b])G 


h+h- 


-b 


—a — b 


—a 


+ exp(— m[a + b])G 


h + h_ 


-2a + 6 


—a + b 


—a 


— exp(— in[a — b])G 


h~h+ 


-b 


—a — b 


—a 


+ exp(— m[a + b])G 


h_h + 


-b 


+a-b 


+a - 2b 


— exp(+i7r[a — b])G 


hrhr 


-b 


—a — b 


—a 


+ exp(— in[a + b])G 


h-h- 


+b 


+a + b 


+a 


— exp(+i7r[a + b])G 





y < z 


z < y < x 


x < y 


I 


h+h+ 


-b 





—a 


-G 


h + h + 


-2a - b 


-2a - 2b 


-a - 2b 


+ exp(-2i7r[a + b])G 


h+h- 


-b 





—a 


-G 


h + h- 


-2a + 6 


-2a 


—a 


+ exp(— 2ma)G 


h-h+ 


-b 





—a 


-G 


h_h + 


-b 


-2b 


+a - 2b 


+ exp(-2i7r6)G 


hrhr 


-b 





—a 


—G 




+b 





+a 


+G 



Table 3: The phases of products h(a,x — y)h(b,y — z) and resulting integral / (last 
column) for the case x < z (upper part) and x > z (lower part) 
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In this way we have proven the composition relation ( |3.12j ). (Note that a = a + 1, b = 
+ 1). 

Alternatively, composition relation can be easily proved in representation given by Eq. 
( ggp for a, (3 < 0. The relations 



/-t-oo 
- - z)dy, a, (3 < (A.17) 

•00 

after inserting from ( |2.28| ) and simple substitutions immediately follow from ( |A.6| ). 
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